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We report results based on a 1.86 fb 1 data sample collected on the T(5S) resonance by the 
Belle detector at the KEKB asymmetric energy e + e~ collider. The inclusive production of 
D s , D° and J/tji mesons at the T(5S) is studied. From the D s inclusive branching fraction the 
ratio f a = (16.4 ± 1.4 ± 4.1)% of B^B^ to the total bb quark pair production at the T(5S) 
energy is obtained in a model dependent way. The exclusive decays B s — > D a *^ + tt~ (/ p~) and 
B s — > J/ty(j)(/ri) are studied and a significant B s signal is observed combining these modes. 
The B s meson production is found to proceed predominantly through the creation of B*B* 
pairs. Upper limits on B s — > K + K~ , B s — > (jyy, B s — » 77 and B s — » D s *^ + D s *^~ decays are 
reported. 



1 Introduction 

The possibility of studying B s decays at very high luminosity e + e~ colliders running at the 
energy of the T(5S) resonance has been discussed in several theoretical papers > 2 . In 2003 
the CLEO experiment collected 0.42 fb -1 at the T(5S) and observed evidence for B s meson 
production in both inclusive 3 and exclusive 4 modes. However, simple calculations assuming an 
approximate SU(3) symmetry indicate that many interesting B s measurements require a data 
sample of at least ~ 20 fb -1 , which is well within the capability of the e + e~ B factories. To test 
the experimental feasibility of such measurements a data sample of 1.86 fb -1 was recently taken 
at the T(5S) energy with the Belle detector 5 at the KEKB 6 asymmetric energy e + e~ collider. 
This data sample is more than four times larger than the CLEO dataset at the T(5S). 

An energy scan was performed just before the T(5S) data taking to define the position of 
the T(5S) resonance production maximum. An integrated luminosity of ~30pb~ 1 was collected 
at five values of energy between 10825 MeV and 10905 MeV, at intervals of 20 MeV. Finally, 
the data sample of 1.86 fb -1 was taken at the T(5S) energy of ~10869 MeV. The experimental 
conditions of data taking at T(5S) are exactly the same as in T(4S) or continuum runs. The 



data sample of 3.67 fb taken in the continuum 60 MeV in CM energy below the T(4S) is used 
to estimate continuum backgrounds. 

2 Number of bb Events 

In the energy region of the T(5S) the hadronic events can be classified into three physics cate- 
gories: mm, dd, ss, cc continuum events, bb continuum events and T(5S) events. The bb continuum 
and the T(5S) events always produce final states with a pair of B or B s mesons and, therefore, 
cannot be topologically separated. We define the bb continuum and T(5S) events collectively as 
bb events. All bb events are expected to pass into the BB, BB*, B*B*, BB-k, BBtttt, BB*tt, 
B S B S , B S B* or B*B* final states. The excited B mesons decay to ground states via B* — > B^y 
and B* — > Bgj decay channels 7 . 

The uu, dd, ss, cc continuum subtraction method is applied to obtain the number of bb events 
in the T(5S) data sample: 



Here N|| is the number of bb events in the T(5S) data sample, and Ngg^ and N^t are 
the numbers of hadronic events in the T(5S) and continuum data samples, respectively. The 
efficiency to select a bb event in the T(5S) data sample, e^g = (99 ± 1)%, and the efficiency ratio 
for hadronic events in the T(5S) and continuum data samples, e5g d V e cont* = 1-007 ± 0.003, are 
obtained from MC simulation. The integrated luminosity ratio ^ss/Axmt is calculated using a 
standard Belle luminosity measurement procedure based on the measurement of Bhabha events. 
The ratio C§s/ C cont = 0.4740 ± 0.0019 is obtained after detailed calculations. From Eq. (1) we 
obtain the number of bb events in the T(5S) data sample, N|| = (5.61 ±0.03 s t a t ±0.29 SJ/S t) x 10 5 . 
The full systematic uncertainty of ~ 5% includes all systematic errors on parameters used in 
Eq. (1) and is dominated by the uncertainty on the luminosity ratio definition. 

3 Inclusive D s Production 

The method of inclusive D s analysis of the T(5S) data sample was developed in 3 . Almost 
the same technique is applied in this analysis. The production of D s mesons at the T(5S) is 
studied in the clean decay mode, Df — > (f)ir + . The D s signals in the T(5S) and continuum 
data samples are seen in Fig. la for the normalized D s momentum region x(D s ) < 0.5, where 
a bb contribution is expected. The normalized momentum is defined as x(h) = P(h)/ P max (h), 
where P(h) and P m ax(h) are the measured momentum and the maximum possible momentum 
of particle h, respectively. Continuum distributions are normalized to T(5S) distributions using 
the energy corrected luminosity ratio. To extract the number of D s events, the candidate mass 
distribution is fitted to a Gaussian to describe the signal and a linear function to describe the 
background. The x(D s ) distributions are shown in Fig. lb for the T(5S) and continuum data 
samples. These two distributions agree well in the region x(D s ) > 0.5, where bb events cannot 
contribute. An excess of events in the region x(D s ) < 0.5 corresponds to inclusive D s production 
in the bb events. 

Subtracting the continuum contribution, applying the efficiency correction, and summing 
over bins within the interval x{D s ) < 0.5, the inclusive branching fraction i3(Y(5S)— > D s X)/2 = 
(22.6 ± 1.2 ± 2.8)% is obtained. This value is significantly larger than the branching fraction 
B(B -> D S X), which has been measured by the CLEO Collaboration 3 to be (9.0 ± 0.3 ± 1.4)% 
and by the BaBar Collaboration 8 to be (8.94±0.16=L 1.12)% (corrected for our value of B(Df — > 
07r + )). The significant increase of D s production at T(5S), compared with T(4S), indicates a 
sizable B s production rate. 




(1) 




M((|>7i + ) (GeV/c 2 ) 



Figure 1: The D s signal (a) at the region x(D 3 ) < 0.5 for the T(5S) (full circles with error bars) and continuum 
(hatched histogram) data samples, and the D s normalized momentum distributions (b) for the T(5S) (full circles 
with error bars) and continuum (hatched histogram) data samples are shown. 



The ratio f s of B s * B s * events to all bb events at the T(5S) can be obtained from equation: 

B(T(5S) -> D s X)/2 = f s ■ B(B S - D S X) + (1 - f s ) ■ B(B - D S X) (2) 

Using B{B S — > D S X) = (92 ±11)% obtained in the reference 3 in a model dependent way 
the value f s = (16.4 ± 1.4 ± 4.1)% is obtained from Eq. (2). This value agrees well with 
f s = (16.0 ± 2.6 ± 5.8)% obtained by the CLEO collaboration 3 from the D s inclusive analysis. 
The dominant contributions to systematic uncertainty are the uncertainty on B(Df — > 4>tt + ), the 
uncertainty on the number of bb events and the uncertainty on the model dependent assumption 
for B{B S -> D,X). 



4 Inclusive D° and J/rp Production 

The inclusive production of D° mesons at the T(5S) is studied in the decay mode D° — ► K~ir + 
(Fig. 2a,b). Large D° signals are seen in Fig. 2a for the T(5S) and continuum data samples 
for the normalized momentum region x(D ) < 0.5. The number of D° mesons as a function 
of the normalized momentum x(D°) is shown in Fig. 2b for the T(5S) and continuum data 
samples. In a manner similar to the D s inclusive analysis, the inclusive branching fraction 
£(T(5SH D°X)/2 = (53.3 ± 2.0 ± 2.9)% is determined. Using the inclusive D° production 
branching fraction of the T(5S), B, and B s decays and replacing D s by D° in Eq.(2), the ratio 



fs 



(18.7 ± 3.6 ± 6.7)% of B^bP events to all bb events at the T(5S) is obtained. This value 
agrees with our D s measurement. 

The inclusive production of J/ip mesons is studied in the decay mode J ftp — > fJ> + fi~ (Fig. 2c, d). 
The J/i/j production at the T(5S) can be estimated in a way similar to that used in the inclusive 
D s analysis. In the region x(J/ip) < 0.5 a prominent J/ip signal (Fig. 2c) is seen for the T(5S) 
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Figure 2: The D° (a) and J/ip (c) signals at the region x < 0.5 for the T(5S) (full circles with error bars) and 
continuum (hatched histogram) data samples, and the D° (b) and J/tp (d) normalized momentum distributions 
for the T(5S) (full circles with error bars) and continuum (hatched histogram) data samples are shown. 



data sample, whereas the J/ip signal in continuum is small. The J/ip normalized momentum 
distributions are shown in Fig. 2d for the T(5S) and continuum data samples. 

Using the PDG value 7 B(J/tp — > /i + n~) = (5.88 ± 0.10)%, the inclusive branching frac- 
tion fi(T(5S)-> J/^X)/2 = (1.068 ± 0.086 ± 0.057)% is obtained by summing events in the 
background subtracted and efficiency corrected J/ip normalized momentum distribution. The 
T(5S) branching fraction can be compared to the B decay branching fraction 7 B(B — > J/tpX) = 
(1.094 ± 0.032)%, because the inclusive J/ip production rate in B and B s decays is expected to 
be approximately equal. 



5 Exclusive B s — > DW+ir", B s — > DW+p" and B s -» J/^f(p(/rj) Decays 



In this analysis we studied the six conventional B s decay modes B s — > Dfw~, B s — > D* + 7r~, 
£? s — > Dfp~, B s — > D* + p~, B s — > and l? s — > J/^rj, which have large reconstruction 

efficiencies and are described by unsuppressed conventional tree diagrams. mesons are re- 
constructed in the </>7r + , and -fLg-ftT + decay channels. To suppress continuum background 
we exploit topological cuts. 

The signals can be observed using two variables: the energy difference AE = E^f — E^ m 

JOM 



and beam-constrained mass Mb c 



momentum of the B s candidate in the e + e center-of-mass (CM) system and E^ m is the 
CM beam energy. The B s mesons can be produced at the T(5S) energy via the intermediate 

with B* — > B s j. Most B s signal events are concentrated within the 



b{*^ Bi*^ channels, 



and B S B S 



ellipsoidal regions, which are kinematically well separated for the B*B^ 
channels. 

The distribution of data in M^ c and AE for the B s — > Dfir~ decay mode is shown in 
Fig. 3a. Nine events are observed within the -B s signal ellipsoidal region corresponding to B*B* 
pair production channel. Background outside the signal regions is small. The Mb c and AE 
scatter plots are also obtained for the B s — > D* + ir~ (Fig. 3b), B s — > Di*^ + p~ (Fig. 3c) and 
-B s — > J /^><j)/{rj) (Fig. 3d) decays. One of the observed S s — > J/^^ candidates is reconstructed 
in the J/* — > ^ + ^ _ mode and one in the J/^ — > e + e~ mode. One candidate is observed in the 
J/^rj decay mode. 
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Figure 3: The M bc and AE scatter plots for the B s -> D+tt (a), B s -> D* + tt (b) and B s -> D { 3 * )+ p (c) and 

B s — + J/^<f)(rf) (d) decay modes are shown. 



The six B s modes shown in Fig. 3 are combined to increase the statistical significance of 
the B s signal. Distributions in AE are obtained separately for events from three Mb c intervals, 
corresponding to B s production proceeding through the B*B*, B*B S + B S B* or B S B S channels, 
respectively. Each of these three distributions is fitted by the sum of a Gaussian to describe 
the signal and a linear function to describe background. The fits yield 20.0 ± 4.8 events and 



1.3 ± 2.0 events for the B*B* and B*B S + B S B* channels, respectively; no events are observed 



in the B S B S channel. 



From these numbers we obtain the ratio a(e + t 



B* s B* s )/ai 



e ' e 



0.94lSno at the T(5S) energy. Potential models 9 ' 10 predict the fraction of B*B* 



r(*) r(*) 

D s D s j — u.ct_ 00 g 

channel over all these channels to be around 70%. 

The B s and B* masses can be extracted from the M^ c fits in the B*B* channel, selecting 
candidates within the -0.08 < AE < -0.02 MeV range. The M bc distribution (Fig. 4a) is fitted 
by the sum of a Gaussian to describe the signal and the so-called ARGUS function 11 to describe 
background. The fit yields the mass value M(B*) = (5418 ± 1) MeV/c 2 . The observed width of 
the B* signal is (3.6 ±0.6) MeV/c 2 and agrees with the value obtained from the MC simulation, 
which assumes zero natural width. 
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Figure 4: The B* (a) and B s (b) mass distributions for events within the —0.08 < AE < —0.02 McV interval, 
corresponding to the B*B* channel. Curves are obtained from the fits described in the text. 



Using events from the B*B* channel we can obtain also the B s mass (Fig. 4b), replacing 
the energy E^J^ m by the term £^ e ^ m - < AE > in the mass formula. The fit shown in Fig. 4b 
yields the B s mass M{B S ) = (5370 ± 1 ± 3)MeV/c 2 and width a(B s ) = (3.6 ± 0.6)MeV/c 2 . 
The second uncertainty in the B s mass value is the systematic uncertainty due to the statistical 
uncertainty on the < AE > measurement. The obtained B s mass agrees with the most recent 
CDF measurement, M(B S ) = (5366.0 ± 0.8) MeV/c 2 . 



6 Rare B s Decays 

The distributions in data of Mb c and AE are obtained for reconstructed B s — > 77 (Fig. 5a), 
B s -» 07 (Fig. 5b), B s -> K+K~ (Fig. 5c) and B s -> D { * )+ Di* ] ~ (Fig. 5d) candidates. Only the 
B s signal regions corresponding to the dominant B*B* channel are indicated in Fig. 5. These 
regions are wider for the B s — > ^7 and B s — > 77 decays, where the energy losses due to photon 
radiation lead to a large tail at lower values of AE. 
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Figure 5: The scatter plots in M bc and AE for the B 3 — > 77 (a), B s — » 07 (b), B s — * K + K (c) and B, 
Di*^ (d) decays. The ellipses are indicating the B a signal regions for the B*B* channel. 



The numbers of events within the signal regions, estimated background contributions, ef- 
ficiencies, and upper limits for these decays are listed in Table 1. For comparison, the PDG 
upper limits are also shown. The number of initial B*B* pairs is obtained by multiplying the 

number of B^ B^ pairs defined in the inclusive analysis by the production ratio of B*B* pairs 
to Bg Bs pairs obtained in analysis of conventional exclusive B s decays. The efficiencies are 
determined from MC simulation. 



Table 1: The number of events in the signal region, estimated background contribution, efficiencies, upper limits 
and PDG upper limits for the B s — + 77, B s — ► cfyy, B s — ► K + K~ and B s — > Di*^~ decay modes. 



Decay mode 


Yield, ev. 


Backg., ev. 


Eff. (%) 


UL (10~ 4 ) 


PDG UL (10~ 4 ) 


B s -»■ 77 





0.5 


20.0 


0.56 


1.48 


B s ^07 


1 


0.15 


5.9 


4.1 


1.2 


B s -► K+K- 


2 


0.14 


9.5 


3.4 


0.59 


B s - D+D- 





0.02 


0.020 


710. 




B s - D*+D~ 


1 


0.01 


0.0099 


1270. 




B s - D*+D*- 





<0.01 


0.0052 


2730. 





The obtained upper limit for the B s — ► 77 decay is about three times smaller, than the 
current best world value. Within the Standard Model this decay is expected to proceed via 
an intrinsic penguin diagram 12 and the branching fraction is expected to be (0.5 — 1.0) x 10 -6 . 
However that branching fraction is sensitive to some Beyond Standard Model (BSM) contribu- 
tions and can be higher by one to two orders of magnitude in some BSM models 13,14 . Important 
results can be obtained using the B s — > decay modes 15 . These modes are expected 

to be predominantly CP eigenstates and, because of expected large branching fractions, should 
lead to a sizable lifetime difference between CP-odd and CP-even B s mesons. Assuming the sum 
of these branching fractions to be around (5-8)%, we can expect about 5 events in each of four 
possible final states with statistics of ~30fb _1 . Within the SM framework such measurement 
can provide an important constraint on the value of Ar s /r s . 
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